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MODELING  CLIMATOLOGY  OP  AREAL  COVERAGE 


Irving  I .  Crlngorten 

Air  Pare*  Oeophyslcs  Laboratory  (LYT) 
ifcnacoa  APB,  HA  01731 


1  .  INTRODUCTION 

Throughout  a  given  area,  at  any  lnatant  In 
tlae,  the  weather,  auch  aa  rainfall  depleted  on  a 
radar  acreen,  or  clouds  In  a  satellite  picture, 
produces  a  pattern  that  Is  unique  to  that  instant 
in  tlae,  not  likely  to  be  duplicated  at  another 
Instant.  The  change,  from  one  Instant  to  anothar, 
is  likewise  unique.  Both  the  pattern  and  lta 
change,  however,  are  governed  by  basic  restraints 
that  characterize  the  cllsate  of  the  area.  The 
relevant  climatology  aust  consist  of  the  atatistlca 
of  the  events  or  the  probabilities  of  events  in 
the  area . 

If  the  probability  Is  known,  say,  of  2A-hour 
rainfall  exceeding  10  aa  at  a  single  observational 
point,  we  aust  ask:  what  Is  the  corresponding 
probability  that  the  saae  quantity  will  be  exceeded 
everywhere  In  a  surrounding  area,  or  aoae  fraction 
of  a  surrounding  area?  It  Is  this  type  of  cllaatol- 
ogy  that  is  the  subject  of  this  paper. 

Satellite  pictures,  recording  as  they  do  the 
areal  coverage  of  cloudiness,  eventually  should 
help  to  provide  the  kind  of  cloud  cllaatology 
that  we  seek .  Radar  has  been  i  valuable  asset  In 
recent  years  because  of  the  panoramic  view  that 
Is  provided  of  precipitation  within  a  circular 
area  surrounding  the  radar  station.  Satellite 
and  radar,  however,  are  recent  tools  .  Generally, 
the  existing  clleatlc  records  contain  single- 
instant  observations  of  weather  at  single-point 
weather  stations.  Recorded  rainfall  aeount  la 
that  at  a  rain  gauge .  lbapcrature  la  that  In  the 
thcraoaeter  shelter,  and  so  on.  The  cllaatology 
of  a  city  usually  la  that  at  an  observational 
point  in  the  city. 

Ulth  respect  to  tlae,  aost  observations  night 
be  considered  aa  Instantaneous,  observed  at  re¬ 
gular  Intervals.  But  there  are  also  tlae-lnte- 
grated  observations,  like  24-hour  totals  of  rain¬ 
fall  .  There  are  observations  of  the  average  tea- 
perature,  and  aaxlaua  and  alnlaua  teaperatures  In 
•  24-hour  period.  All  such  observations  contri¬ 
bute  to  the  cllaatologlcal  character  of  a  loca¬ 
tion  as  it  relates  to  tlae  or  the  duration  of 
events,  la  there  a  slallar  dlaatlc  record  of 
•patlal  coverage  of  dlaateT 

In  aost  previous  studies,  correlation  coeff 1- 
tents  are  the  principal  object  of  iavestlgatloa, 

*lth  notable  exceptions.  Schreiner  and  Riedel 
(W7|)  have  collected  rainfall  data  iatagratad 


over  areas  ranging  In  size  froa  26  ka2  to  2600 
k.2,  and  have  found  frequencies  of  extreaes  of 
rainfall  aa  a  function  of  size  of  the  ground  cov¬ 
ered.  Other  studies,  by  Court  (1961),  Briggs 
(1972)  and  Roberta  (1971)  have  linked  the  rain  In 
a  circular  area  to  the  central  slngle-atatlon 
amounts  by  Idealized  models.  In  general,  however, 
efforts  et  areal  coverage  have  been  few  and  aodels 
of  areal  coverage  have  been  rare .  The  following 
exaaples  Illustrate  the  potential  of  such  studies . 


Fig.  1.  A  nap  showing  the  location  of  14$  sta¬ 
tions  equipped  with  rain  gauges  In  New 
England . 


Fig.  1  shows  southern  New  England,  centered 
In  Massachusetts.  There  are  14$  little  circles 
showing  locations  of  stations  that  report  24-hour 
rainfall,  lhere  are  $  nesting  squares  of  areas 
ranging  froa  1000  k a2  to  2d, 000  k*2 .  Some  20 
years  of  record  have  provided  the  estimates  (Table 
l)  of  the  aaxlaua  24-hour  rainfall  amounts  aa  a 
function  of  areal  size.  The  frequency  of  no  rain 
at  a  single  station  la  622,  but  within  a  sizeable 
area  this  frequency  aust  drop  steadily  with  in¬ 
creasing  slse  of  the  srea.  The  slngle-pnlnt  prob¬ 
ability  of  2$  aa  of  rain  is  141,  or  the  probabil¬ 
ity  that  it  will  not  be  exceeded  is  97.22.  As  we 
look  at  areas  of  Increasing  slsa,  this  cumulative 
probability  aust  decrease.  Or  coaversely  the  prab- 
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New  Engisnd  July  maximum  24-hour  rein 
fell.  In  square  areas  ranging  (n  *ize 
iron  1,000  km2  to  28,000  km2. 


igle-  Frequency  cf 

station  Areal  .laximuc 

Frequ.  (Area  In  units  of  102  in2) 


ability  that  rainfall  will  exceed  23  no,  somewhere 
in  an  area,  oust  increase  with  increasing  size  of 
the  area . ' 


-J;^7 


* 

a  /> 


Illustration  of  a  PPI-rader  picture,  in 
wV.ich  there  were  scattered  storas  produc¬ 
ing  radar  echoes  . 


rainfall  echoes.  For  example,  at  Key  West,  Flor¬ 
ida,  in  summer,  the  single-point  probability  of  a 
radar  echo  is  IX  .  But  the  radar  sees  rain  over  an 
area.  We  had  divided  the  radar  picture  into 
cells,  each  cell  of  size  1380  km2.  TWo  cells 
have  twice  that  area,  all  64  cells  have  an  area 


0/(0)  aeeai.  cove* 


AREA  -A-  km 

Fig.  2.  The  cumulative  frequency  of  Che  maximum 
24-hour  rainfall  in  the  areas  of  10*  km2, 
10*  km2  and  2.8x10*  km2,  July  in  New 
England.  The  curves  were  fitted  by  a 
model  of  areal  coverage. 

The  numbers  that  appear  in  Tkble  l  also 
apoaar  in  Fig.  2  at  the  X's.  The  solid  curves 
have  been  drawn  by  means  of  a  modal  of  areal 
cover,  which  I  shall  try  to  describe  shortly.  In 
order  to  fit  the  curves  to  the  observational  data, 
a  value  had  to  be  found  for  the  model's  parameter, 
called  scale  distance  (2.6)  km  in  the  example). 

To  understand  this  parameter,  for  the  present,  we 
describe  it  as  the  distance  over  which  correlation 
coefficient  la  0.99. 

Another  question  lei  What  la  the  probability 
•f  rain  ever  a  fraction  of  the  whole  area?  Pig.  ) 
la  n  eketth  of  one  picture  of  a  radar  PPl-acope, 
ahowftag  eca tiered  arena  of  atony  weather,  fig.  6 
la  a  graphical  preaeatatlea  of  frequent  lea  of 
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Pig.  6.  The  frequencies  of  radar  echoes  covering 
at  least,  (3/10)ths  of  the  areas  ranging 
from  1.6x10*  km*  (l  cell)  to  10s  km2 
(64  cells).  The  plots  are  for  Key  West, 
Florida  and  Cape  Ihtteraa,  N.C.  In 
tuoner  (all  hours  Included). 


of  approximately  100,000  km2.  The  frequency  of 
echoes,  la  at  leant  S/lOths  of  small  aroma,  la 
.  aaarly  the  same  IX.  If  the  area  becomes  large, 
then  a  21  alaglapolat  probability  will  mat 
plaid  SOX  caver  very  often.  Am  solid  curves  at 
thin  diagram,  again,  ara  provided  bp  a  medal  mf 
areal  coverage. 
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Pig.  5.  The  cumulative  probability  of  cloud 

cover,  fro*  clear  to  overcaat,  for  the 
area  of  the  aky  dome  that  an  obeerver 
aeea,  and  for  25  times  that  area.  The 
exaaple  la  for  Geraan  station*  with  scan 
cloud  cover  0.76  In  January. 


fbrtpartlal  cover  let  ua  look  at  fig.  5,  for  - 
an  exaaple  of  cloudiness.  In  this  dlagraa  we  have 
plotted  probabilities  of  cloud  cover  f roa  clear  to 
overcast.  Germany  has  a  satisfying  high  frequency 
of  clouds  In  winter,  averaging  762 .  Du  suaaarles 
of  ground  obaervatlons  at  Bitburg  are  shown  by  the 
X's.  The  model  provided  a  reasonable  fit  to  ’his 
data,  but  the  model  does  more.  Consider  a  larger 
area  that  a  satellite  alght  see.  Obviously  the 
probability  of  all-clear  auat  be  less,  likewise 
of  complete  overcast.  The  larger  area  must  expe¬ 
rience  more  cases  of  partial  cover.  Ihe  aodel 
provide*  auch  estimate*  that  are  otherwise  not 
evaiiable. 

One  characteristic  that  we  must  simulate  with 
°ur  oodel  1*  the  change  of  correlation  coefficient 
with  the  distance  between  stations.  Fig.  6  is  an 
example  taken  from  a  1962  report  (Bortonl  and 
Iwod) .  Correlation  begins  high  for  short  dis¬ 
tance*,  but  drops  to  sero  and  becomes  negative  at 
•ubstantlal  distances  between  stations.  Ihe  sx- 
•■Ple  shows  ce  of  surface  temperatures  plotted 
•gainst  distances  between  stations  from  a  few  km 
to  several  thousand  km  apart. 

.  Viaoratlcally  an  abundance  of  data,  at  an 
•hmi«c«  of  stations  In  a  email  area,  with  a 
*****  euecession  of  observatlona,  could  provide 
the  kind  of  areal  climatology  that  we  seek.  Prac- 

1 colly  such  abundance  would  confront  us  with  the 
••ossoue  task  of  mapping  and  recording  events  as  a 
watery  of  the  weather,  prior  to  Its  suaarlsatloe 


Fig.  6.  Correlation  coefficient  versua  distance 
between  stations,  of  wintertime  surface 
temperatures  (Kef:  Bertonl  and  brad, 
1964). 


In  statistical  format .  this  Is  the  difficulty 
that  we  wish  to  overcome  by  modeling.  He  want* to 
be  able  to  simulate  a  sequence  of  changes  In  the 
weather  covering  an  area,  whenever  we  choose,  and. 
we  want,  ultimately,  to  provide  usable  statlatlc* 
of  areal  coverage.  Including  cc'a,  and  probabili¬ 
ties  of  fractional  cover  by  weather  conditions 
such  as  cloud  cover. 

By  and  large  we  must  depend  on  statistical 
models  which  I  must  yet  describe.  But  there  Is  a 
concept  that  needs  to  be  examined  first,  concern¬ 
ing  the  equivalent  normal  deviate. 


Fig.  7.  Illustration  of  t ran* format Ion  of  a  wea¬ 
ther  element  Into  its  equivalent  normal 
deviate  (or  vice  versa),  (sample  la  for 
January  noontime  visibility  at  Bedford, 
MB. 
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2.  *'  }UIVAL£HT  NVtMAL  DEVIATE  (END) 

Fig.  7  1*  a  **aple  plot  of  the  cumulative 
frequency  of  a  weather  element,  in  this  case  of 
visibility  at  Bedford,  MA,  in  January,  at  noon- 
tioe.  The  X's  mark  the  frequencies  that  have  been 
summarized  froa  20  years  of  data.  They  show,  for 
example,  a  frequency  of  2X  for  visibility  less 
than  1/4  alia,  7X  for  visibility  less  than  1  mile, 
also  62X  frequency  of  visibility  greater  than  10 
miles.  Alongside  the  scale  of  cumulative  prob¬ 
ability  we  can  plot  another  scale  of  the  equiva¬ 
lent  normal  deviate,  (y).  The  correspondence  of 
y  to  F(y)  can  be  found  in  nearly  every  textbook  on 
statistics.  By  means  of  a  diagram  such  as  Fig.  7, 
or  alternatively  by  equations,  we  obtain  a  trans¬ 
formation  of  the  variable  (V  in  this  case)  into 
its  END  (y),  or  vice  versa. 

Our  models  are  derived  using  the  END'S.  The 
results  are  then  transformed  to  the  appropriate 
weather  variables,  when* so  required. 

3.  APPROACH  ID  AREAL  CLIMATOLOGY 

An  effective  model  of  areal  cover  will  pro¬ 
duce  a  simulated  horlaontal  field,  stochastically, 
without  recourse  to  physical  laws  or  dynamics. 

TVio  early  models  (Crlngorten,  1979)  were  devised 
in  the  manner  illustrated  in  Fig.  8.  Imagine 
that  our  horizontal  space  is  filled  with  a  fine- 
mesh  grid  of  points.  Assign,  initially,  random 
numbers  to  all  points.  Around  point  (l,j)  draw  a 
circle  and  take  a  weighted  average  of  all  the 
random  numbers  within  the  circle,  then  replace 
the  random  number  at  (l,j)  with  the  weighted  aver- 
age,  y(J,j).  If  we  do  this  for  all  points  within 
a  large  square,  we  shall  obtain  a  pattern  wherein 
the  new  values,  at  all  points,  will  be  correlated 
with  each  other.  Models  that  have  been  composed 
in  this  way  differ  froa  each  other  in  the  weight¬ 
ing  function  that  is  selected  for  each  model. 

Once  a  model  la  selected,  the  edopted  proce¬ 
dure  can  be  repeated  many  times  to  produce  many 
simulated  synoptic  situations,  thus  comprising  a 
large  sample  that  can  be  surveyed  for  spatial  cor¬ 
relations  or  associations. 

This  early  modeling  procedure  has  been  set 
aside  in  favor  of  a  new  model  that  requires  the, 
generation  of  only  a  few  dozen  random  numbers  per 
map.  Instead  of  many  thousands. 

4.  BOEHI’S  SAWTOOTH  WAVE  MODEL  (BSW) 

The  presently  preferred  model  was  devised  by 
Major  Albert  R.  Boehm,  USAF.  lh  has  given  it  the 
descriptive  name  Sawtooth  Wave.  We  shall  refer 
to  it  as  the  BSW-model . 

We  begin  with  e  horlsontel  spece  (Fig.  9) 
with  axes  (U,V)  and  point  of  origin  (0,0).  We  In¬ 
troduce  a  stationery  wave  formation  with  waves 
chat  are  straight  and  uniform  In  wavelength.  The 
wavelength  (A)  becomes  the  parameter  of  the  model, 
dilek,  la  our  climatic  fields,  measures  several 
hundred  to  several  thousand  km.  Xn  the  following 
derivation,  however,  the  wavelength  is  the  unit 
of  distance. 


Fig.  8.  To  illustrate  the  former  development  of 
stochastic  modeling. 


Fig.  9.  To  illustrate  the  stochastic  procedure 
la  the  BSW-model. 


In  cross-section  the  waves  are  saw-tooth  in 
form.  The  left-hand  edge  of  each  wave  has  zero 
height,  the  right-hand  edge  has  a  height  of  unity. 
At  any  point  (u,v)  in  between  the  wave  phaee, 
x(u,v),  is  uniformly  distributed  between  tero  and 
one. 

The  wavs  formation  hat  two  degreeo  of  freedom, 
end  thus  allows  for  selection  of  two  random  nutr* 
berss  one  for  the  direction  (I)  of  the  wove  for¬ 
malism  and  one  for  the  phono  me  the  point  of  ori¬ 
gin,  n(0,0),  which  moot  ho  ogmml  to  the  fractional 
distance  t  shown  in  the  diagram.  At  ovary  point 
in  the  B-f  field  chore  la  a  wave  phaee,  x(u,v) . 
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.Soy,  consider  aany  (S)  wave  formations,  each 
with  4  random  direction  (0t)  and  each  with  a 
randoc  phase  (>  j )  *t  the  origin.  Then,  at 
each  point  (u,v)  tn-*re  will  be  u  values  of 

*l(u,v)  for  1  «  1  «  .« 

0  «  Oj  <  2» 

0  «  T,  <  1 

The  average,  y(u,v),  of  all  a  values  of  a^(u,v), 
as  N  becoaes  large,  will  have  an  asymptotically 
Gaussian  distribution.  When,  for  convenience,  H 
is  Bade  equal  to  12,  the  equation  becoaes 


y(u,v)  -  1  xg(u.v)  -  6 

There  will  be  a  value  of  y(u,v)  at  each  point 
(u,v)  in  the  horizontal  space.  A  sample  of  such 
generated  fields  Is  shown  in  Fig.  10,  with 
distinct  low  and  high  cells  and  saddle  points. 


Fig.  11.  Craphlcal  presentation  of  the  change 
of  cc  with  distance  In  the  BSU  model . 


4b.  ‘Areal  coverage 

While  we  have  an  analytical  solution  for  the 
correlation  decrease  with  distance  we  lack  such 
good  fortune,  at  thla  time,  for  the  prlnary  goal 
of  this  paper:  to  find  the  probability  of  the  maxl- 
nua  condition  In  an  area  or  a  fraction  of  the  area. 
We  have  obtained  approximate,  although  good,  solu¬ 
tions  by  Monte  Carlo  simulation.  The  following  Is 
a  brief  description  of  our  procedures . 

In  a  field  generated  by  the  model,  consider  a 
square  area  of  side  s  (In  units  of  the  wavelength)  . 
At  some  point  In  the  area  (o2)  there  will  be  a 
maximum.  Them  will  be,  also,  a  tenth  largest,  a 
two-tenths  largest,  —  and  finally  the  minimum  In 
the  area.  In  terms  of  a  fraction  (f)  of  the  area 
the  maximised  minimum  In  such  fraction,  y(f,a2), 
becomes  the  primary  subject.  We  want  to  find  the 
cumulative  probability  distribution  of  the  thres¬ 
hold  value:  F(y:f,a2). 

A  large  number  (M)  of  stochastically  gener¬ 
ated  maps  will  provide  a  frequency  distribution  of 
y(f,s2)  that  will  approach  the  true  probability 
distribution  asymptotically  as  M  Is  Increased. 


Fig.  10.  Example  of  horizontal  field  generated 
stochastically  by  the  BSW-model. 


4a .  Correlation 

An  analytical  solution  has  been  found  for  the 
correlation  coefficient  p(s)  between  points  sepa- 
rated  by  distance  (s),  in  units  of  the  wavelength: 


This  kind  of  Monte  Carlo  simulation  was  con¬ 
ducted  with  If  “  25,000.  On  each  map  12  square 
areas  were  surveyed  for  the  largest  value  of  y, 
the  (1/I0)th  largest,  (2/10)th  largest,  and  so  on, 
down  to  the  smallest  value.  The  linear  size  of 
each  square,  in  units  of  wavelength,  was 

s  ■  2*/340  where  s  ■  -1(1)10 


b(s)  »  3s2  -  (12/*)s  ♦  1  for  0  <  s  <  1 
-  3s2  -  (12/s )s  ♦  1 

♦  (24/e)lcos-l(l/s)  -  / s2— 1 ) 

for  1  <  s  <  2 


From  all  25,000  maps  cumulative  frequencies  were 
found  for  each 

y(f,a2)  -  -3 .5(0. 5)3 .5 


The  graph  of  p(s)  versus  s  is  shown  in  Fig.  11. 
For  our  purposes  we  shall  rarely  need  to  consider 
distances  of  more  than  one  wavelength.  The  dis¬ 
tance  (s)  will  usually  be  a  number  between  zero 
■nd  one,  given  by 

s  •  e*/A 

•'  is  the  distance  (km)  and  A  is  ska  para- 
•atar  wavelength  (km) . 


for  f  -  0(0 .1)1.0 
and  s  -  -1(1)10. 

Doing  the  frequencies  thus  obtained  the  charts  of 
Fig.  12  ware  drawn  to  a  how  the  probability  esti¬ 
mates  of  y-maxlma  for  f  «  0.0,  0.1,  0.3,  0.5 
ooccesalvely,  amd  for  a  from  -1  to  0.  Smoothing 
woe  aaeoaaary,  especially  at  the  ext  rones  of  the 
dlatrlhmtloaa . 


Ac .  Applications 

Por  a  square  area,  A  ka*,  with  linear  dioen- 
alon  s’  ka  we  sat 

a  -  e'/A 

for  which 

s  “  in(340s)/tn  2 

When  a  la  entered  In  a  chare  of  H|.  12,  the 'ver¬ 
tical  line  upward  froa  the  abscissa  provides  eatl- 
aatea  of  the  probabilities  of  various  thresholds 
of  y  corresponding  to  the  area  (A) . 

Instead  of  the  wavelength  (A)  it  Is  conven¬ 
ient  to  substitute  the  paraacter  (r),  called  the 
scale  distance.  In  the  BSW-model  It  is  a  distance 
over  which  the  cc  la  0.9888  (close  to  0.99).  By 
and  large,  it  la  a  aore  aeaningful  paraaetar  In 
application,  given  by 

r  -  s' /340s  ka  -  A/340  ka  -  /A/2*  ka. 

Or,  for  z: 

z  -  tn(/A/r)/tn  2 

For  the  New  England  24-hour  wlntar  precipitation, 
we  found  scale  distances  of  9  to  10  ka.  for  radar 
echoes  they  ranged  froa  1  to  2  ka  In  suaaer,  2  1/2 
to  4  1/2  ka  in  winter.  Cloud  cover  scale  dis¬ 
tances  have  ranged  froa  1/2  ka  to  10  ka  depending 
on  season,  tlae  of  day  as  well  as  geographic  loca¬ 
tion.  * 

5.  TIME  CHANGES 

The  above-described  aodel  enables  us  to  gen¬ 
erate  'snapshot”  pictures  stochastically  of  a 
horizontal  field .  So,  far,  the  ENl>-value  y(u,v) 
at  a  point  (u,v)  on  one  aap  is  Independent  of 
that  generated  on  any  other  aap. 

Let  us  suppose  that  one  aap  has  been  gener¬ 
ated  at  tlae  t  with  values  of  yt(u,v)  at  each 
point  (u,v)  of  the  space,  let  another  aap  be 
generated,  for  the  tlae  (t+4t),  In  terns  of 
an  END  (n)  to  give  values  of  n(t+gt\(u,v) 
at  each  point  (u,v).  Make  the  codification  to 
y(t+6t)(u,v)  by  aeans  of  the  Markov  aodel*. 

5r(t+«t)("',,)"‘>t  •J't(u,',,‘H^l~pt  ’  n(t+«t)<“*T* 
where  the  cc  (pt)  for  tlae  lag  (<t)  is  given  by 

Pt  •  ozpf —  t/T  ) 

"here  T  (In  hours)  Is  a  paraaetar  called  relaxa¬ 
tion  tlae. 

Between  the  y-values  of  the  later  aap,  the 
co  can  be  proved  to  be  still  p(a),  the  aaae  as 
It  was  on  the  aap  at  tlae  t.  that  la,  the  cc- 
tield  an  the  new  aap.  It  hours  later,  la  pre- 
Mr*sd.  Vila  la  a  necessary  condition  to  assure 
continuity  of  the  nodal's  eharactorlatlca .  An 
**^fle  of  hourly  changes  la  shown  In  figs.  1)  for 
*  -  20  hours. 


Fig.  13.  Ezaaple  of  stochastic  horizontal  fields 
following  one  another  at  24-hour  Inter¬ 
vals  In  a  Markov  process . 


6.  REMAINING  PROBLEMS 

These  efforts  have  raised  aore  probleas  than 
they  have  solved.  Me  have  devised  a ‘procedure  to 
generate  synoptic  pictures  of  the  weather,  stochas¬ 
tically,  without  regard  for  the  underlying  phys¬ 
ical  processes.  He  have  aade  no  effort  to  use  the 
aodels  of  Nuaerlcal  Heather  Prediction  (NHP)  to 
produce  successive  pictures  of  the  aeteorologlcal 
events.  At  this  tlae,  however,  we  do  not  know  If 
the  hourly  sequences  of  MMP  will,  in  fact,  yield 
a  cllaatology  that  approxlaates  that  of  tha  real 
dlaate  with  raspect  to  such  elcaents  as  cloud 
cover,  especially  aesoscale  features.  At  this 
point  In  tlae  we  are  justified  In  adopting  a  pro¬ 
cedure  that  will  produce  aaps  that,  collectively, 
will  yield  a  cllaatology  consisting  of  a  good 
approxlaatlon  to  the  horizontal  field  of  correla¬ 
tion  and  satisfactory  probability  estlaatcs  of 
the  aaslaua  and  alnlaua  of  the  weather  conditions 
In  fraction  (f)  of  an  area  (A).  There  la. nothing 
In  the  aodel,  yet,  that  will  cause  a  synoptic  sit¬ 
uation  to  be  advected  as  well  ss  nod if led. 

An  ultlaate  aodel  should  be  4-dlaenelonal,  to 
yield  probsbllltiss  of  conditions  In  3-diaen- 
slonal  space  and  to  Include  eventa  In  tine.  More¬ 
over,  tbs  probeblllty  of  events  or  conditions 
should  be  estlaated  analytically  instaad  of  by 
Monte  Carlo  s inula t ion.  Ms  have  aade  a  beglan- 

but  ouch  raaalns  to  be  done. 
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